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ABSTRACT 



Or 



Context. The optical counterpart of the neutron star X-ray binary and well known X-ray burster, 4U 1636-53 (= 4U 1636-536 = 
V801 Ara) has been well studied for three decades. However to date, no infrared studies have been reported. 
Aims. Our aims are to identify and investigate the near-infrared (NIR) counterpart of 4U 1636-53. 
Methods. We present deep, AVband (2.2yum) imaging of the region of 4U 1636-53 taken with the Infrared Spectrometer And Array 
Camera (ISAAC) on the Very Large Telescope. Archival optical and UV data are used to infer the 0.2 - 2.2yum spectral energy 
■ distribution (SED). 

Results. One star is located at coordinates a = 16:40:55.57, 6 = -53:45:05.2 (J2000; lcr positional uncertainty of ~ 0.3 arcsec) 
which is consistent with the known optical position of 4U 1636-53; its magnitude is K$ = 16.14 ± 0.12. This star is also detected 
in the 2MASS survey in /-band and has a magnitude of J = 16.65 ± 0.22. Under the assumption that the persistent emission is 
largely unvarying, the 0.4 - 2.2yum de-reddened SED can be described by a power law; F v oc v ' 5±03 , with some possible curvature 
(F v oc v sl 5 ) at 0.2 - 0.4/im. The SED can be approximated by a blackbody of temperature ~ 27000 K. This is typical for an active 
OS ', low -mass X-ray binary, and the emission can be explained by the outer regions of a (likely irradiated) accretion disc. We therefore 

C^*) . interpret this As -band star as the NIR counterpart. 
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1. Introduction radiated accr etion discs which also dominate at optical wave- 
lengths (e.g.. Ivan Paradijs & McChntockll 1995b . 

Low mass X-ray binaries (LMXBs) are interacting binaries 
where a low mass donor is transferring material onto a neu- 
tron star or a black hole. In order to transport the excess of Even in the near-infrared (NIR), where the companion star 
angular momentum outwards, an accretion disc is formed. In could have an important contribution, the disc emission seems 
' the disc, the gravitational potential energy is transformed into dominant in bright systems. For instance, in the prototypical per- 
mainly X-ray radiation and kinetic energy, and temperatures sistent neutron star, Sco X-l, which given its relatively long 
approach ~ 10 8 K. The mass transfer rate supplied by the orbital period (~ 19 h) should have a large, evolved compan- 
donor star, Mi, is driven by the binary /donor evolution and ion star, no spectral feature from the donor has been detected in 
mass transfer rates Mi > M CI it ~ 1O~ 9 M kg yr~' (where M the NIR (Bandvopadhva v et al.1 1997b . In most persistent neutron 
is the mass of the Sun in kg) result in persistently bright X- star X-ray binaries, spectral and temporal studies have favoured 
ray sources dKing et alJ Il996l) . There are ~ 200 such bright an X-ray heated accr etion disc as the or igin of the NIR emission 
(L x - 10 36 - 10 38 erg s" 1 ) LMXBs in the Galaxy and most of (e.g., 4U 1705-440; THoman e t al. 2009). Compact jets produc- 
them harbour neutron stars as implied by the detection of pul- ing synchrotron emission typically dominate the radio emission 
sations and X-ray bursts resulting from nuclear burning caused (Migl iari & Fenderl l2006) and their spectra can extend to higher 
by the accumulation of Hydrogen and Helium on their surfaces, frequencies. In the N IR, high amplitude flares from GX 17+2 
They show energy spectra dominated by emission from their ir- (Call anan et alJl2002l) . a synchrotron spectrum in 4U 0614+09 

(Migliari et al.ll2010h and v ariable linear polarization in Sco X- 

* Based on observations collected at the European Southern 1 (Russell & Fender 2008) have suggested a strong infrared jet 

Observatory, Chile, under ESO Programme ID 085.D-0456(D). contribution in these persistent neutron star X-ray binaries. 
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Fig. 1. VLT / ISAAC deep, high resolution /T s -band finding chart 
for 4U 1636-53. The optical position (USNO-B1 0.339" error 
circle) is indicated by a white circle and the infrared (ISAAC) 
position is shown by two black lines. The two positions agree 
within lcr errors. The six 2MASS stars are indicated by diagonal 
lines. 

One of the classical neutron star systems is 4U 1636-53 (= 
4U 1636-536 = V801 Ara). It is an X-ray burster, which has 
been extensively studied in X-ray and optical regimes for more 
than three d ecades (e.g., Pedersen et al. l 1982). It has a 3.8 h or- 
bital period (van Paradii s"et al.lll990HGiles et al.ll2002h pointing 
to a relatively faint, late type companion star. The spectrum of 
the system from X-ray to optical wavelengths seems to be domi- 
nated by the emission of its bright accretion disc, the companion 
star only being detected by using emission lines arising from re- 
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processing of the strong X-ray emission (L x 
in its inner hemisphere (Casares et al. 2006). 

To date, no detections of 4U 1636-53 have been reported at 
wavelengths longer than the optical regime. At radio frequen- 
cies, upper limits of both the pe rsistent and burst fluxes were 
presented in lThomas et alj dl979l) . Here we present the first de- 
tection of the near-infrared counterpart of 4U 1636-53. Together 
with available optical and UV data, we construct the intrinsic 
0.2-2.2 //m spectral energy distribution (SED). 

2. Very Large Telescope observations 

The data were acquired with the Infrared Spectrometer And 
Array Camera (ISAAC) on the European Southern Observatory 
(ESO) 8-m class Very Large Telescope UT3 (Melipal) on 2010- 
05-21 02:23 - 03:17 UT (MJD 55337.118 + 0.019) under ESO 
Programme ID 085.D-0456(D). Twenty-four data cubes, each of 
exposure time 132.5 s (comprising 530 x 0.25 s individual ex- 
posures) were obtained of the region of 4U 1636-53 in high 
time resolution (FastJitter) mode using the A^-band filter cen- 
tred at 2.16yum. The brightest object in the field has < 1000 
counts in each 0.25 s exposure, so non-linearities are not a con- 
cern. A small (up to 6") telescope offset was applied between 
cubes. Average images were made of each cube, removing the 
sky generated from a median of the surrounding averaged cubes. 
The field of view of each ISAAC image is 40" x 40" and the 



pixel scale is 0.1478" pixel . The 24 average images were then 
aligned and stacked using IRAF to produce a deep image. The 
total on-source exposure time of the combined image is 3180 s. 
The As -band seeing as measured from the combined image was 
0.6 arcsec and conditions were clear. The camera was rotated at 
an angle of 127° to maximise the number of reference stars on 
the array. The combined image is presented in Fig.Q]and can be 
used as a deep, high resolution A"s-band finding chart. 



3. Source identification 

Within the combined imag e, six stars listed in th e Two Micron 
All Sky Survey (2MASS; ISkrutskie etal]|2006l) are detected, 
which were used to calculate the world coordinate system and 
achieve £"s-band flux calibration. DAOPHOT in IRAF is used to 
perform point-spread-function (PSF) photometry on the com- 
bined image. All comparison stars were successfully fitted by a 
single PSF consistent with the image PSF; no blends were iden- 
tified. There is one star consistent with the best optical position 
of 4U 1636-53, as listed in the USNO-B1 digital sky survey 
{a = 16:40:55.61, 5 = -53:45:05.1; J2000, with an error circle 
of 0.339"), shown as a white circle in Fig.[T] We derive a position 
of (x — 16:40:55.57, 6 = -53:45:05.2 (J2000) for the infrared 
counterpart of 4U 1636-53, with a lcr positional uncertainty of 
~ 0.3 ar csec. Other optical coordinates of 4U 1636-53 in the lit- 
erature dJernigan et al.ll 1 977b Eiuet al.ll200ll:ISamus et all 2003) 
are also consistent with the USNO-B1 and ISAAC positions. No 
precise X-ray coordinates of 4U 1636-53 from Chandra or Swift 
observations exist in the literature. 

Some faint, blended stars lie just 1.5" from 4U 1636-53 in 
the ISAAC image, but none are consistent with the optical posi- 
tion. The above source is the only candidate counterpart in our 
T^s-band image. All stars within a 3" radius of the LMXB are 
> 1 .2 magnitudes fainter in /fs-band than this proposed counter- 
part. No faint blended stars are visible after subtracting the PSF 
of the counterpart. 

From the known magnitudes of the 2MASS stars within the 
field, we measure the magnitude of 4U 1636-53 during the time 
of our observations to be Ks = 16.14 + 0.12. The errors are 
dominated by systematic errors derived from the six 2MASS star 
magnitudes. 4U 1636-53 is detected with a signal-to-noise ratio 
(S/N) of 145. The limiting magnitude in the image is Ks ~ 19.3 
(i.e. stars brighter than this are detected at the 3<x level). 

The counterpart we have found is not listed in the 2MASS 
point source catalogue. On inspection of the 2MASS images, 
the counterpart is not visible in K and //-bands, but appears as 
a faint source in /-band. Performing PSF photometry on the J- 
band 2MASS field, we obtain a magnitude of J = 16.65 + 0.22 
for 4U 1636-53; the S/N is 7. The aforementioned faint stars 
within ~ 2" of the source may contribute up to ~ 10 - 30% of 
the flux measured, although the PSF fitting method should have 
removed the majority of this excess flux (which is offset from 
the central PSF position). This 2MASS observation was made on 
1999-06-18 04:27 UT. In the 2MASS //-band image, the upper 
limit to the magnitude of 4U 1636-53 is H > 15.76. The new 
VLT and 2MASS detections of 4U 1636-53 are listed in Table 
1. 



4. Swift UVOT observations 

Optical-ultraviolet observations of 4U 1636 -53 were made with 
the Ultra Violet/Optical Telescope (UVOT; iRoming et al]|2005l) 
on board the Swift satellite and the data are publicly available. 
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Fig. 2. Intrinsic (de-reddened) UV-optical-NIR spectral energy 
distribution of 4U 1636-53 (top) and observed (reddened) fluxes 
(lower, crosses). The systematic errors due to uncertainties in the 
extinction are demonstrated by the dotted error bars. 

Table 1. New NIR - optical - UV magnitudes of 4U 1636-53. 
The Swift UVOT photometry is taken between 2005-02-13 and 
2010-06-10. 



Filter X 1 S/N - Magnitude Days 2 





(A) 




mean 


brightest 


faintest 






21590 


145 


16.14 


16.14(12) 


16.14(12) 




J 


12350 


7 


16.65 


16.65(22) 


16.65(22) 




V 


5468 


28 


17.94 


17.94(6) 


17.94(6) 


1 


b 


4392 


21 


18.67 


18.67(8) 


18.67(8) 


1 


u 


3465 


27^3 


17.82 


17.63(7) 


17.95(8) 


5 


uwl 


2600 


10-24 


18.35 


18.16(9) 


18.67(15) 


4 


um2 


2246 


7-15 


18.95 


18.61(13) 


19.31(20) 


5 


uw2 


1928 


7-34 


19.20 


18.53(10) 


19.61(21) 
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'The central wavelength of the filter in Angstroms; 2 The number of 
dates the source was observed with UVOT. 



UVOT observed 4U 1636-53 on 11 dates between 2005-02-13 
and 2010-06-10. On most dates, one or two of the six UVOT fil- 
ters were used. The image data of each filter on each date were 
summed using uvotimsum. Photometry of the source in indi- 
vidual sequences was derived with uvotsource. 4U 1636-53 
is clearly detected at the position derived above, in all six fil- 
ters (5500 - 1900A). An extraction region of radius 3" centred 



at the ISAAC-derived position was adopted, which excludes a 
neighbouring star 6" away from 4U 1636-53; no stars appear to 
contaminate the flux of 4U 1636-53 within the extraction radius. 
The mean magnitudes and range in each filter are given in Table 
1. 



5. Intrinsic spectral energy distribution 

We de-redden the /-band and UVOT magnitudes using the 
known interstellar extinction towa rds 4U 1636-53; Ay = 2.5 + 
0.3 mag (derived from o ptical data;lLawrence et alj!983l) apply- 
ing the extinction law of ICardelli et all & 19891) . We also take op- 
tical m agnitudes from Lawre nce et alj ( 1 1983b and Pedersen et alj 
(fT982h : U = 17.92 + 0.06, B = 18.47 + 0.03, V = 17.89 + 0.03, 
R = 17.46 + 0.03, / = 16.77 + 0.03, and de-redden them 
in the same manner. Using the de-reddened fluxes (adopting 
the standard conversion zeropoints for optical Johnson filters, 
and 2MASS IR filter s from the Explanationary Supplement of 
Skru tskie et al.l 120061) we construct the intrinsic NIR-optical- 
UV SED of the persistent (non-burst) emission from 4U 1636- 
53. 

The SED, which spans one order of magnitude in frequency, 
is presented in Fig. [2] The observed, reddened fluxes are plotted 
in addition to the de-redde n ed dat a. Th e optical (/ to 1 7-band ) 
SED from lLawrence etail (|1983l) and iPedersen etal] d 19821) . 
which has small errors, has a spectral index of a = 1.5 + 0.4, 
where F v cc V . The error on a is dominated by the uncertainty 
in the extinction. This spectral index is typical for active low- 
mass X-ray binaries (both neutron st ar and black hole systems; 
e.g.. lHvnesll2005l : lRussell et alJl2007h and is consisten t with the 
outer regions of a blue accretion disc. lShih et al.l ( 1201 lb showed 
that the optical emission is correlated with the soft X-ray flux; 
the optical light is dominated by the irradiated disc (X-ray repro- 
cessing on the disc surface). 

The K§, //-band and /-band magnitudes correspond to de- 
reddened flux densities of F Y ^ S - 0.30 + 0.03 mJy, F Vy n < 0.83 
mJy and F v j = 0.67 + 0.14 mJy, respectively. The optical- 
infrared (Ks to //-band) SED (neglecting the UVOT data) can 
be fitted with a power law with spectral index a — 1.5 + 0.3 
(solid line in Fig. [2J. This is similar to the optical spectral in- 
dex, and the Z^-band and /-band fluxes are consistent with an 
extrapolation of the ac cretion disc spectrum. It is u nlikely that 
IR emission lines (e.g. Bandvopadhva vet alJll997h could con- 
tribute a significant fraction of the observed IR fluxes. There is 
some possible curvature in the spectrum at the higher frequen- 
cies, apparent in the UVOT data, and the whole SED can be 
approximated by a blackbody at a temperature of ~ 27 000 K 
(dotted line in Fig. [2]). The curvature implies the blackbody of 
the irradiated disc spectrum may peak at around 1900 A, sim- 
ilar to some other LMXBs dHvnesI 12005). However, the curva- 
ture (and even more so the spectral index of the power law) 
is sensitive to the uncertainty in the extinction, and the optical 
counterpart varies by ~ 0.6 mag amplitud e (a factor of ~ 1.7 
in flux) over weeks to months (e.g.. IShih et al.l [20 1 lb . We de- 
tect variability in our UVOT data over several years (Table 1). 
Quasi-simultaneous optical-infrared data are required to accu- 
rately measure the spectral index, and uncertainties will be de- 
creased once the extinction towards the source is measured more 
accurately. 

There is no evidence for the companion star, or synchrotron 
emission from the jet (if the system has a jet) to dominate the 
infrared flux. Both of these components are expected to be red- 
der than the irradiated disc component, and would produce an 
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infrared excess above the disc spectrum. Since the spectrum be- van Paradijs, J. & McClintock, J. E. 1995, in Optical and Ultraviolet 
tween J and Xs-bands is blue (a ~ 1.5), these components make Observations of X-ray Binaries, eds, Lewin, W. H. G., van Paradijs, J., van 

a negligible contribution to /-band. In the most extreme see- den Heuvel - E - R J - (Cambridge, Cambridge Univ. Press), 58 
nario, the disc could have a steep spectrum, a — 2.0 between 
J and TTs-band, and therefore the disc could be as faint as 0.18 
mJy in /Ts-band (0. 12 + 0.03 mJy fainter than the observed flux). 
If this were the case, the jet or companion could contribute up to 
46% of the A^-band flux. These calculations again neglect pos- 
sible variability, since the NIR observations were made on dif- 
ferent dates. Simultaneous optical-infrared observations could 
constrain the level of star or jet emission more precisely. 

6. Conclusions 

We have identified the near-infrared counterpart of the neutron 
star burster LMXB 4U 1636-53. Its magnitudes on the dates ob- 
served are K s = 16.14 + 0.12 (VLT / ISAAC), J = 16.65 + 0.22 
(2MASS). The intrinsic infrared-optic al-UV spectrum of the 
persistent (non-burst) emission is consistent with a blackbody, 
likely from the irradiated surface of the accretion disc. We find 
no evidence for emission from other components that may be ex- 
pected to contribute, such as the donor star or synchrotron emis- 
sion from a jet, although simultaneous optical and infrared data 
are needed to constrain these contributions further. 
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